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Abstract 

We study the spectrum of the reggeized gluon compound states in high-energy QCD. 
The energy of the states depends on the total set of quantum numbers whose values are 
constrained by the quantization conditions. To establish their explicit form we apply the 
methods of nonlinear WKB analysis to the Schrodinger equation for the N reggeon state 
in the separated coordinates. We solve the quantization conditions for the N = 3 reggeon 
states in the leading order of the WKB expansion and observe a good agreement of the 
obtained spectrum of quantum numbers with available exact solutions. 
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1. Introduction 



In perturbative QCD analysis of the Regge asymptotics of hadronic scattering amplitudes, the 
power rise with energy of the physical cross-sections is attributed to the contribution of the color- 
singlet compound gluonic states. These states are built from an arbitrary conserved number, 
N — 2, 3, of reggeized gluons, or reggeons, and they can be defined as solutions of the effective 
(2 + 1)— dimensional Schrodinger equation [l[ 

7-L N \XN,{q}) = E Ny { q y\xN,{q}} ■ (1-1) 

The effective QCD Hamiltonian T~Cn describes a pair- wise interaction between N reggeons on the 
2-dimensional plane of their transverse coordinates, p& = (p Xk ,p yk ) with k = 1, N, and its 
explicit form can be found in M . The longitudinal gluonic degrees of freedom define the rapidity 



of reggeons that enter into (|1.1| ) as a evolution time. The wave function XN,{q}(pi ~Po, Pn~ po) 
depends on some set of quantum numbers {q} and the center-of-mass coordinate po of the reggeon 
compound state. 

The contribution of the iV reggeon states to the physical cross section has a form 



cr is) ~ s 



and at large energies s it is dominated by the state with the maximal energy, E = max^} Ex,{ q }- 
Therefore solving the Schrodinger equation ( |1 . 1|) one is mostly interesting in finding the state 
with the energy Eq, or equivalently the ground state of the hamiltonian —Hn- 

The Schrodinger equation ( |1 . 1|) has a number of remarkable properties. The effective QCD 
Hamiltonian is invariant under projective SX(2,C) transformations of the 2-dimensional 
gluon coordinates pj 



with ad — he = 1 and (zj, Zj) being (anti)holomorphic reggeon coordinates. Denoting the gener- 
ators of these transformations as 



N N N 



^z k d k , L_ = ^(9 fc , L+ = ^z 2 k d k} (1.3) 



k=l k=l k=l 



and having similar expressions for the antiholomorphic generators, one may verify that [Hn, L a ] = 
[Hn, L a ] = 0. This symmetry allows to classify possible solutions to (|1 . 1| ) according to the 
representations of the <SX(2,C) group. In particular, the wave function is transformed under 
projective transformations as 

X N ,{ q} -> (cz + d) 2h (cz + df h XN){q} (1.4) 

with (zq,zo) being (anti)holomorphic coordinates of the center-of-mass. The conformal weights 
of the state, h and h, take the values corresponding to the unitary principal series representation 
of the SL(2,€) group § 

1 ~h Tl — 1 — Tl 

h = h = Yiv, {n = 7L, v = R) (1.5) 
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parameterized by the pair (n, v) of integer n and real v. These parameters determine the spin 
s = n and the scaling dimension d = l + 2iv of the state and control the transformation properties 
of the state under rotations, z k — > z k , z k — > e _l * z k , xn —> e m<?i 'Xn an d dilatations, z k — > \z k , 
z k — > Xz k , xn —> A ~ 2w Xn- 

Being expressed in terms of holomorphic and antiholomorphic reggeon coordinates, the effec- 
tive Hamiltonian splits into the sum of two mutually commuting (1 + 1)— dimensional Hamilto- 
nians acting separately on (anti)holomorphic coordinates, Hn = Hn + Hn and [Hn,Hn] = 
[j|]. This property is exact for iV = 2 and N = 3 reggeon states and for iV > 4 it holds in 
the multi-color limit of QCD. It implies that the original equation can be replaced by the 
system of two (1 + 1)— dimensional Schrodinger equations in holomorphic and antiholomorphic 
subspaces. As a result, the wave function of the state can be factorized as 



Xn{z,z) 



22C a ip N}0 ,(z)(p N ,g(z) , (1.6) 



where z = {z\ — Zq, zn — zq) denotes the set of A" holomorphic coordinates and the coefficients 
C a are fixed by the condition for Xn(z, z) to be a well defined function of its arguments on the 
2-dimensional plane of transverse coordinates. Additionally, the energy En of the state can be 
evaluated as the sum of holomorphic and antiholomorphic energies that one finds solving two 
(1 + 1)— dimensional Schrodinger equations. 

It turns out that the holomorphic Hamiltonian Hn coincides with the Hamiltonian of 1- 
dimensional Heisenberg magnet of noncompact SX(2,(D) spin zero p. There is however an 
important difference between the holomorphic Hamiltonian Hn and well-known integrable (1 + 
1)— dimensional Hamiltonians 0,0- Examining the properties of Hn one finds [5] that in contrast 
with the full (2 + 1)— dimensional effective Hamiltonian, the operator Hn is unbounded on the 
space of the wave functions ( |1.4j ). The same property holds for the antiholomorphic hamiltonian 
Hn- It is only the sum of holomorphic and antiholomorphic (1 + 1)— dimensional Hamiltonians 
that is bounded from above. This happens because the holomorphic and antiholomorphic energies 
depend on the common set of quantum numbers ( |1.5|) and infinities cancel in their sum. 

Analyzing the Schrodinger equation (|1.1|) one finds || || that it possesses the family of 
mutually commuting holomorphic and antiholomorphic integrals of motion, {q k , q k }, k = 2, N 
defined as 

q k = i k Yl z hh- z 3kh d h- d 3k > [q k ,H N } = 0, (1.7) 

• * ./ •••• jr -V 

where zj k = Zj — z k and the definition of the antiholomorphic charges q k is similar with Zj replaced 
by Zj. The total number of the conserved charges, 2(A^ — 1), is large enough for the A^— gluon 



Schrodinger equation (|1 . 1| ) to be completely integrable. The wave function of the A^— reggeon 
state, XN,{ q }, diagonalizes these charges and one identifies the corresponding eigenvalues as a 
complete set of quantum numbers {q} introduced before in (|i~T|). One also takes into account 
that on the space of the wave functions (|1.4j) the holomorphic and antiholomorphic charges are 
related as || 

q k = q* k , k = 2,...,N. 

The integrability of the Schrodinger equation ( |1 . 1|) ensures that the effective Hamiltonian H,n 
acting on the 2A^ coordinates of reggeons can be expressed as a complicated function of 2(A^ — 1) 
integrals of motion (q k , q k ) plus 2 operators reserved for the components of the center-of-mass 
momentum of the state. We notice however that (anti) holomorphic components of the total 
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momentum coincide with the generators, L_ and L_, of the projective transformations ( |1.3|) and 
the dependence of the Hamiltonian on them is protected by the SL(2,€) symmetry, [Hn, L a ] = 0. 
This leads to the following identity [[| 

En,{q} = £jv(<?2, <?3, •••) Qn) + £jvfe, #3, Qn) , 

where En denotes the energy of the (anti)holomorphic Hamiltonian. 

Thus, the problem of calculating the spectrum of the effective QCD Hamiltonian (|1.1|) is 
reduced, first, to finding the explicit form of e n as a real function of the conserved charges q 2 , 
qiv and, second, to establishing the quantization conditions for their eigenvalues. The first part 



of the problem was studied using different approaches [|T0|, |IT] and similar results for were 
obtained. In the present paper we address the problem of finding the spectrum of the conserved 
charges. 

The paper is organized as follows. In the Sect. 2 we employ the WKB expansion to obtain 
the quantization conditions for the eigenvalues of conserved charges. We show in Sect. 3 that 
these conditions have the form of the Whitham equations. Solving the Whitham equations in 
the special case of N = 3 reggeon states we obtain the quantized charges q 2 and q-$ and study 
their properties in Sect. 4. Concluding remarks can be found in Sect. 5. Appendix contains some 
properties of the elliptic integrals which are helpful in solving the Whitham equation. 



2. Quantization conditions 

The "lowest" integrals of motion, q 2 and q 2 , provide the simplest example of the quantization 
conditions. One verifies using ( |1.3| ) and ( |1.7| ) that q 2 and q 2 coincide with the quadratic Casimir 
operators of the SL(2,<C) group, q 2 = Lq(Lq + 1) + L^L + , and for the wave functions transforming 
as ( |1.4j ) their eigenvalues take the form 

q2 = -h(h-l) = --[~ + iv) , q 2 = -h(h - 1) = q* 2 , (2.1) 
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being discrete in n and continuous in v. The quantization conditions for "higher" integrals of 
motion q^ and qt (k > 3) are more involved |], [T(J. One can get some insight into their properties 



by exploring the relation between the wave functions of reggeon compound states xn(z,z) and 
2-dimensional conformal field theories. To this end, let us consider the following function || |l(| 



Xn{z,z) = (z 12 z 23 ...z N1 ) 1 xn{z,z) 



According to its properties under the projective transformations, (|1.4|), the function xn can be 
represented as (N + 1)— point correlation function in 2-dimensional conformal field theory 

Xn(z,z) = ((f)(z 1 ,z 1 )...(j)(z N ,z N )O h - h (z ,z )) . (2.2) 

Here, (p(z, z) is a quasiprimary reggeon field with the conformal weights (1, 0) and O h f l (zo, z ) is 
a quasiprimary operator interpolating the reggeon compound state. The operator O h % has the 
conformal weights ( |1.5| ) and it is transformed under ( |1.2| ) according to the irreducible principal 
series representation (irreps) of the SL(2,€) group labeled as t^ 2n \ For the reggeon operator 
(j)(z,z) the corresponding irreps can be identified as t^°' 2 \ The product of iV reggeon operators 
in (|2.2|) is reducible with respect to the action of the SX(2,(D) transformations. It corresponds 
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to the tensor product of N copies of t*- ' 2 -* which in turn can be decomposed into irreducible 
components using the following fusion rules || 



t (0,2) ^(0,2) = Q t (u,2n) t t (v,2n)® t m= f( 



f.(p,2m) 



In conformal field theory these relations are implemented in the operator product expansion of 
the quasiprimary operators (f>(z, z) and O h j l {z, z) ]12], 0. Applying the fusion rules N — 1 times 
one obtains an expansion of the product (f>(zi, zi)...(f>(zN, %n) in (^3) into an infinite sum of the 
irreducible components. Each component can be parameterized in two equivalent ways: either 
by the conformal weights (h, h) plus the set of 2(N — 2) quantum numbers (qk, qk), or by the set 
of 2(iV — 1) numbers (n^, Vk) corresponding to the iV — 1 intermediate states t^ Uk ' 2nk \ This leads 
to the following dependence 



qk = 5fc((ni, vi), (njv-i,^jv-i)) , (k = 2, N) (2.3) 



together with h = 1+n £- 1 s r i VN _ x ari d h — 1 — h*. Thus, apart from the set (tojv-i, "n-i) defining 
the conformal weight of the N reggeon state, one expects the appearance of (N — 2) additional 
sets (rife, Vk) parameterizing the eigenvalues of the operators q% and q~k- In particular, for given 
conformal weights (h, h) the possible values of the quantum numbers qk and qk form the family 
of continuous functions of real u\, vn-2, different members of which are labeled by integers 

Let us consider two special cases. In the first case, one takes all v\~ = and rik = odd. 
Using (|1.5| ) one finds that it corresponds to restricting the conformal weights of the intermediate 
states t^fe ,2n fe) and the iV reggeon state, t ( - UN - 1 ' 2nN ~ 1 \ to take only integer values. The spectrum 
of the operators qk and qk is expected to be discrete. Indeed, this result was previously obtained 
within the algebraic Bethe Ansatz approach |J by numerical solution of the Bethe equations || 
and later the analytical expressions for quantized qk and qk were derived in the quasiclassical 
approximation [0. In the second case, one chooses nk = and Vk are arbitrary. The conformal 
weights of the intermediate states have the form hk = \ + ivk and they run along imaginary 
axis. It has been argued |§ that it is for these values of the conformal weights that the energy of 
the N reggeon states approaches its maximal value. The quantized qk and qk become continuous 



functions of real Ui, Un-i- As was shown in |14|| , these functions can be found within the WKB 



approach as solutions of the Whitham equations. In the next Section we derive these equations, 
solve them and study the properties of quantized qk and g&. 



2.1. Separated variables 

The eigenvalue problem for the operators qk and q~k defined in (|1.7| ) leads to a complicated system 
of partial differential equations of the N— th order on the wave function xn- One can significantly 
simplify the problem by using the fact that quantization of the charges qk and qk does not depend 
on the particular representation in which the wave function xn is defined in order to go from 
(z, z)— coordinates to a new set of separated holomorphic reggeon coordinates 0, ||, xi,...xn, and 
their antiholomorphic counterparts Xk = x* k . 

In the separated coordinates the wave function of the N reggeon state is factorized into the 
product of functions each depending on only one coordinate 

Xn(x, x) = Q(xi)...Q(x N -i) x N h x Q(xi)...Q(x N -i) x^ h . (2.4) 
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Here, the function Q(x) is a solution of the second order finite-difference equation || 

Q(x + i) +Q(x-i) = x~ N A(x)Q{x), A(x) = 2x N + q 2 x N ~ 2 + ... + q N (2.5) 

known as the Baxter equation and Q(x) satisfies similar relation with all replaced by % = q%. 

The separated gluon coordinates have the following properties. The coordinates xn and xn 
correspond to the center-of-mass motion and they are defined as eigenvalues of the operators 

xn = LqLZ 1 , xn = LqLz 1 (2.6) 



with L_ and L being the SX(2,<D) generators (|1.3|). The definition of the remaining coordinates 



Xk and Xk for k = 1, .., iV — 1 is more involved and it can be found in [/J. It will be important 
for us that they satisfy the relations |7j 

[L ,x k ] = [L_,x k ] = 0, [L , x n ] = ix N , [L-,x N ] = -i 

with k — 1, N — 1 and similar relations hold for Xk and xn- From these relations one finds 
that the separated coordinates Xk and x^, being functions of the gluon (zj, %)— coordinates, are 
invariant under translations, — > Z\~ + e, and dilatations, z& — > Xz k , while the coordinates xn 
and xn are transformed in the same way as z k and Zk- Therefore, in the expression (|2.4j ) for the 
wave function in separated coordinates, it is the factor x N h x N h that carries the spin n = h — h and 
the scaling dimension d — 1 + 2iv — h + h of the reggeon state. Moreover, making the shift, 
Zk — ► -Zfc — zq, one can restore the dependence of the wave function (|2.4j ) on the center-of-mass 
coordinate of the state, (z , z ), as xn — > xn — z and Xk — > x k (k — 1, N — 1). 

The wave function ( p..6|) has to be a well defined function of the reggeon coordinates z k and 
Zk- Performing transition to the separated coordinates one imposes the same condition on ( |2.4Q 
as a function of the complex x k and = x^.. Let us first examine the dependence of \n on 
the coordinate xn- If Xn is a single valued function of complex xn it should not be changed 
as Xn ^ e 2m XN in ( |2.4|) . One checks that this is true provided that the spin h — h = n 
of the state is integer. Requiring additionally that the scaling dimension of the state has to 
bed = h + h = l + 2iv one reproduces the quantization conditions (|1.5|) for the conformal 



weights. We notice that although the wave function (p.4j) stays invariant under 2ir rotations on 
the Xn— plane, its holomorphic and antiholomorphic components get nontrivial monodromies for 
all values of the conformal weights h and h except for integer ones. The latter case corresponds 



to the polynomial solutions studied in 0, ITU 



To satisfy the condition of singlevaluedness for the dependence of ( [2.4|) on the remaining 
coordinates Xk one has to examine the properties of the product Q(x)Q(x*) on the complex 
x— plane. As we will show below, this consideration leads to the selection rules for the quantum 
numbers g 3 , qN 

2.2. Properties of the Baxter equation 



The function Q(x) entering ( [2.4|) obeys the Baxter equation (|2.5| ). This equation has two linear 
independent solutions, Q+(x) and Q-(x), satisfying the normalization condition 

W(x) = Q + (x + i)Q-(x) -Q-(x + i)Q+{x) = 1 

and having the following asymptotics at infinity 

Q+{x) ~ x h , Q-{x)~x l ~ h asx^oo. (2.7) 
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The analysis of the Baxter equation for the function Q(x) is similar since the functions (Q(x))* 
and Q(x*) satisfy the same equation. As a result, one gets the following general expression 

Q(x)Q(x*) = CapQa(x)(Qp(x))* 

a,P=± 

with C a/ 3 being arbitrary coefficients. However, it follows from ( |2.7|) that for the conformal 
weights h and h of the form ( J1.5| ) the functions Q+(x) and Q-(x) acquire nontrivial monodromies 
at infinity as x — > e 2m x, which should cancel in the product Q(x)Q(x*). This requirement leads 
to C = C ++ = and 

Q(x)Q(x*) = C + -Q+(x)(Q-(x)Y + C. + Q.(x)(Q+(x)y . (2.8) 

with C + _ and C_ + being some constants. 

The asymptotic behaviour ( |2.7| ) of the functions Q±(x) is fixed by the conformal weight of 
the state and it is not sensitive to the value of "higher" quantum numbers q^, q^. To obtain 
the quantization conditions for q 3 , q^ one has to identify the remaining singularities of the 
solutions of the Baxter equation, Q+(x) and Q-(x), on the complex x— plane and examine the 
monodromy of Q+(x)(Q-(x))* around them. 

As an example, let us seek for solutions Q{x) to the Baxter equation within the class of 
rational functions. The distinguished property of such solutions is that they have a finite number 
of poles on the x— plane and scale at infinity as an integer power of x. According to (|2.7| ) the 
latter property implies integer values of the conformal weight h. Then, straightforward analysis 
shows that Q(x) can not have poles on the x— plane, since otherwise in order to balance the both 
sides of the Baxter equation ( J2.5| ) the functions Q(x) should have an infinite sequence of poles 
shifted by ±i along imaginary axis. In similar manner, taking the limit x — > in ( |2.5| ) one can 
argue that Q(x) should vanish at the origin faster then x N . These two conditions fix Q(x) to be 
a polynomial of degree h > N with N— times degenerate root x = 

h-N 

Q+(x) = x N l[(x - A fc ) . (2.9) 

k=l 

The remaining h — N roots of the polynomial, satisfy the Bethe equations and for given 
h they take a discrete set of real values || Substituting this solution back into the Baxter 
equation Q2.5p and comparing its both sides one can express quantized qk in terms of the roots. 
The spectrum of qk takes the form (|2.3|) with Vf. = and n& = odd and it can be calculated 
explicitly ||. Thus defined polynomial solutions of the Baxter equation were first obtained within 
the algebraic Bethe Ansatz approach and their properties were studied in detail ||, p~0|] . 

To find the most general form of the quantization conditions ( |2.3|) one has to solve the Baxter 
equation ( |2.5|) for arbitrary complex conformal weights of the form flL5| ). Different attempts 
have been undertaken |15, [L0[ 0, but the analytical solution similar to ( p.9|) is not available 



at present. As the first approximation to the exact solution, it has been proposed jnj to look 
for the asymptotic solution to the Baxter equation by applying the methods of nonlinear WKB 
analysis well-known from the studies of the Painleve type-I equation |T7 . 



2.3. WKB expansion 

Let us interpret the Baxter equation ( |2.5| ) as a 1-dimensional discretized Schrodinger equation 
on the complex x— line and let us look for its solutions in the form of the WKB expansion. 
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One observes that the Planck constant is formally unity K = 1 in the Baxter equation and the 
quasiclassical limit corresponds to large values of the "energies" q k entering into the definition 
( p.5| ) of A(x). This limit can be effectively performed by introducing an arbitrary small parameter 
r) and rescaling the coordinate, x — > x/rj, and the quantum numbers, q k — > qkV k -> m (HID- Then, 



the asymptotic solution to the Baxter equation has the form |T^, [U], [Hfl 

Q(x/rj)=exp(^J dS(x)j , S(x) = S (x) + ijS^x) + 0{rf) , (2.10) 
where S (x) and Si(x) satisfy the relations 

2cosh^(x) = ^, S[(x) = i -^-\nsmhS' (x) (2.11) 

CC Zi doc 

with S'(x) = d x S(x) and 

A(x) = 2x N + q 2 x N ~ 2 + ... + q N , q k = q k 7] k = fixed as r? -> . (2.12) 

Similar expressions can be obtained for the antiholomorphic solution Q(x). 

The limit 77 — > 0, or equivalently = large, is an analog of the classical limit in quantum 
mechanics. As 77 — » the quantum fluctuations become frozen and the motion of particles is 
restricted to their classical trajectories defined by the leading term in (|2.10f) , the eikonal phase 



Sq(x). Having identified So(xk) as the "action" function for the system of N reggeons in the 
separated coordinates one may find the classical trajectories of gluons generated by the "classical" 
Hamiltonians q k by solving the system of the corresponding Hamilton- Jacobi equations. This 
system of equations is exactly integrable and its solutions can be described as cycles on the 



Riemann surface Tn in the following way ]TS, O 
Let us introduce the algebraic curve Tn 



V: y 2 = A 2 (x)-Ax 2N , (2.13) 



with A(x) given by Q2.12 ). This relation defines y = y(x) as a double- valued function on the 



complex x— plane having 2A^ — 2 branching points at x — <jj such that y(<Jj) = or 

l 2 (a 3 )=Ao 2N . (2.14) 

The same function becomes a single- valued on the Riemann surface Tn obtained by gluing 
together two sheets of the x— plane along the cuts foj, &2j+i], J = 1, AT — 1 running between 
the branching points. Encircling the cuts by oriented closed contours one constructs the set of 
N — 1 cycles a%, ckjv-i- The genus of T N equals N — 2 and it depends on the number of 
reggeons "inside" the compound state. 

The classical trajectories of reggeons in the separated coordinates x\, xn-i correspond 
to the independent motion along A^ — 1 different a— cycles on the Riemann surface T^- The 
branching points of become the turning points of the reggeon trajectories. Coming back from 
the separated coordinates to the original reggeon coordinates, and one can show that the 
same motion looks on the 2-dimensional plane of transverse gluon coordinates as a soliton wave 
propagating in the system of A^ interacting gluons [|14]]. Its explicit expression was constructed 



in [2(| using the methods of the finite-gap theory out of the curve Tn. The quantum numbers 
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of the N reggeon state, q k , become parameters of the soliton solutions. The quantization of q k 
appears as a result of imposing the Bohr-Sommerfeld quantization conditions on the periodic 
reggeon trajectories in the action- angle variables (a*., (p k ) defined as |L9| 



1 

2^ 



dS (x) 



dSp(x) 
da k 



(2.15) 



where < (p k < 2ir and a k are cycles on the Riemann surface Q2.13| ). 

To establish the explicit form of the WKB quantization conditions let us consider the prop- 
erties of the differential dS = dS (x) + rjdSi(x) + 0{rf) defined by the WKB ansatz ( |2.10| ). One 
finds from ( |2.11 ) that dSo is the following meromorphic differential on 



dSo(x) = dxlnu(x) 



dx 

y 



NA(x) - xA'(x) 



(2.16) 



Here the notation was introduced for the curve uj = u(x) 



uj + 1/u = A{x)/x 



N 



x N (u 



l/u) 



(2.17) 



and the sign ~ means equivalence up to the exact differential, d(x\nu>), having vanishing 
a— periods. 

Solving ( 2.11]) we obtain two different solutions, dSo t +(x) and dSo,-(x), which are continu- 
ous functions of complex x except across the cuts [<J2ji 02j+i], satisfying the relation dSo, + (x) = 
—dSo-(x). Being combined together they define the differential dSo(x) as a double- valued func- 
tion on the complex x— plane. Nevertheless, dSo is a single-valued function on the Riemann 
surface built from two copies of the a;— plane to which we will refer as upper and lower sheets. 
Examining (|2.11 ) as x —>■ oo we obtain the behaviour of dSo(x) on the upper, dSo i+ , and lower, 



dSo _, sheets of T N as 



dS ±(x) 



x 



oo 



(2.18) 



where q2 is given by Q2.ll) . 

The first nonleading correction to the WKB phase S(x) is defined by the differential dS\ 
which one calculates using ( |2.16|) and (|2.17|) as 



i A\x)-Ax 2n 
dSi{x) = ~d\n ^ 



2N-2 




x - a k 



(2.19) 



with o"fc being the branching points ( j2.14|) . In contrast with dSo, this differential is well defined 



on the complex a;— plane. It has poles located at the branching points (|2.14j ), at the origin and 
at the infinity where its asymptotics is given by 



dSi(x) 

2 x 



as x 



oo . 



(2.20) 



Combining together ( p. 16 ) and ( 2.19| ) we find that two different branches of the differential dSo 
give rise to two linear independent WKB solutions to the Baxter equation 



Q±(x/rj) = exp 



dS±(x] 
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where 



dS + (x) = dS 0i+ (x) + r)dSx(x) + 0(r] 2 ) , 

dS-(x) = dS Q ,-(x) + V dSi(x) + 0(rj 2 ) = -dS 0t+ (x) + TjdSt + 0(r?) 



Their asymptotics at infinity can be obtained from ( 2.18|) and ( |2.20|) as 



dx 



ds ± ( x ) = -iii— (^±v^2 + - + o(n) 

We recall that at small rj the quantum numbers q k scale as q^ ~ l/i] k . Therefore, replacing 
y/—q2 = h — | + 0(h~ 2 ) and taking the limit 77 — ^ one checks that the functions Q±{x/r]) 
behave at infinity in accordance with ( |2.7|) . 

Substituting the obtained WKB expressions for Q±{x/rj) into ( |2.8| ) and keeping only the first 
term in the r.h.s. of ( |2.8|) one finds the holomorphic and antiholomorphic WKB solutions to the 
Baxter equation as 

Q{x/rj) = ex p{^y dS ,+ + iJ dS 1 + 0{ n )j (2.21) 

Q(x*/r}) = expj^y^^ -i^dS^j +0(rj) 

Let us examine the analytical properties of these functions on the complex x— plane. Relations 
ffngD and (l27T9f) imply that the singularities of both functions are located on the cuts [cr 2 k, C2fc+i]- 



For the product Q{x/rj)Q(x* /rj) to be a uniform function of x one has to require that for x 
encircling the cuts [<7 2 fc, c^+i] 011 the complex plane along the cycles it should come to the 
starting point with the same value. This gives the following set of iV — 1 quantization conditions 

Re ^ dS ,+J - *7 Im (j> dS i^j + °(V 2 ) = V™ k (2.22) 
where n k (k = 1, ...,N — 1) are integers and 0(rf) stands for nonleading terms in the WKB 



expansion. According to ( 2.19 ), the differential dS\ has poles on the x— plane and the calculation 
of its a— periods is reduced to taking the residues at the poles located inside a k . Since by 
definition each a— cycle encircles two branching points, o 2 k and cr 2 k+i, and, additionally, one of 
the cycles, say ak , encircles the pole at x — 0, one gets 

If 1 N 

-j dS 1 = -- + 6 UM -. (2.23) 



Substituting this expression into ( 2.22 ) we write the quantization condition in the following form 



i £ dSo,+{x) =v(y + Wk ) + ' ( 2 - 24 ) 

with Pf. (k = 1, N — 1) being arbitrary real numbers. This relation has a simple interpre- 
tation in terms of the action variables ( 2.15Q . Namely, it states that the possible values of the 



action variables a k match into the spectrum of the conformal weights ( |1 . 5| ) of the principal series 
representation of the SL(2,G) group. 
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The periods of the differential dSo,+ entering the l.h.s. of the quantization conditions (|2.24j) 



are uniquely defined by the values of the quantum numbers q2, qN- Therefore, solving the 
system of N — 1 equations ( |2.24| ), one will get q 2 , qN as functions of N — 1 complex numbers 
— ^2 + ^ V hi each corresponding to N — 1 cycles on the Riemann surface Tn- This result is in 
agreement with our expectations ( |2.3| ) and, additionally, it implies that the WKB quantized q k 
are holomorphic functions of l k 

qk = Qk{ti,—,tN-i), 4=y+^fc- (2.25) 

The same parameters determine the value of the conformal weight h as follows. One takes into 
account that the sum of the a— periods of the differential dS 0i+ is equal to its period around the 
infinity where the differential has a simple pole ( |2.18|) with the residue t]\/—q2- In this way one 



gets using fl2.24j ) 



JV-l 



(t + iVk ) + 0(,,) (2 ' 26) 

k=l 

with r] ~ l/y/—q2. For large — q2, or equivalently large conformal weights h, one replaces a/ - <?2 = 
h — 5 + 0(1/ h) to obtain the expression for h which reproduces possible quantized values of the 
conformal weight ( |TT5| ) with n = J2k n k an d v = J2k v k- 

Having defined in ( |2.23| ) and (|2.24j) the periods of the differentials, one calculates the mon- 



odromy of the WKB solutions to the Baxter equation, Q2.21|) , along the a^— cycles as 

g (x) g (x) eM^+^+T**.*.) , Q(x*) x ^ Q(x*) e M^-^~^o) , ( 2 .27) 

where higher order WKB corrections were neglected. One verifies that the monodromies cancel 
against each other in the product Q(x)Q(x*). Moreover, moving x around all a— cycles on the 
complex x— plane and subsequently applying ( p.27|) one reproduces the monodromy at infinity 

Q(x) -> Q(x) e 2mh , Q(x*) -> Q(x*) e 2 ^ l ~ h *) . 

Relations (|2.27|) allow us to establish the correspondence between WKB quantization conditions 
( |2.24| ) and the special class of exact polynomial solutions to the Baxter equation ( |2.9| ). Namely, 
comparing ( |2.27|) with (pT9 ) one has to require that e^p(2m(^^- + iv^ + Y^k,k )) = 1- This leads 



to the additional constraints on the parameters of the WKB quantization conditions ( |2.24 ) 

Vk = , n k + 5 kM N = 2Z+ + 1 
and, as a consequence, the conformal weight h = | + ^ . ^ + ivj takes only positive integer 



values. The resulting quantization conditions were studied in |14| and their solutions, q k , 
were found to be in a good agreement with numerical calculations within the algebraic Bethe 
Ansatz approach M. 



3. Whitham flow 



Let us apply the quantization conditions (|2.24j) to obtain the expressions for the charges qs, q^ 
in the leading order of the WKB expansion. In the leading order one neglects O (rj 2 )— corrections 
to (pi) and uses invariance of the differential (|2.16|) and the curve ( |2.13|) under transformations 



x 



Xx . 



qk ^ Xq k 



dSo — > XdSo 
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induced by reparameterization of the WKB expansion parameter, 77 — > Xrj, to put rj = l/^/—q 2 , 
or equivalently q\ = — 1, in (|2.24j) . This amounts to replacing q\ in the definition, fl2.12|) and 



( |2.16| ), of the differential dSo by invariant parameters 

u k = qk/(-Q2) k/2 , k = 3,...,N, (3.1) 



called the moduli of the curve T^. One also introduces notation for the r.h.s. of (|2.24j ) as 



^ + j Vl N - 1 

2 E** = 1 ' ( 3 - 2 ) 



k=l 

where the last relation follows from (12.261). 



In terms of new variables the quantization conditions (|2.24D become the relations between 



the moduli Uk and the parameters 5j. In particular, in the case of N = 3 states the quantization 
conditions for the moduli 113 = u look like 



(3u — 2x)dx 

2tt4, k — 1,2. (3.3) 



Qfc a/(u — x) (4a; 3 — x + u) 



Solving (|3.3j ) and (|2.24|) one finds the moduli, Uk = Uk(5i, ...,Sn-i), as smooth functions of the 



complex parameters 5k and then calculates the corresponding values of the quantum numbers q%, 
q N using ( |OD and (|3~^ ). 

Let us consider the relations (|3.1|) and fl3.2|) in two special cases discussed before Sect. 2.1: 
rife = 0, z/fc 7^ and rife 7^ 0, = for all k, to which we will refer below as I and II, respectively. 
We find that in both cases the parameters (|3.2|) take real values 



A (i) _ Vk _ n k 



Ui + ... + v N _i n 1 + ... + n N _i 

They are discrete for and continuous for 8% with the former being the subset of the latter. 
Calculating the moduli as u k = Uk{5^) and u k = Uk{5^) one realizes that they enjoy the 
same properties. Among other things this means that in both cases the moduli are universal, 
that is they are given by the values of same function u = u(5) on the real 5— axis. In particular, 
this function can be calculated for real 5 in the case II using available exact polynomial solutions 
of the Baxter equation and then applied to evaluate the quantized qk in the case I. 

Universality of moduli leads to the following expressions for the quantum numbers ( |3.1| ) and 

m 

$ = g« = M^,...,^) (3.4) 



and 



with v = J2k=i u k and n = Ylk=i n k- One notices that the substitution 



maps one of the relations into another one and the same transformation leaves invariant the 
expression (|2.25|) and the quantization conditions ( |2.24j ). 
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We would like to stress that the relations (|3.4j) and (|3.5|) were obtained in the leading order 
of the WKB expansion and they are valid for large n and v. In this limit, one finds using (|2.1|) 
that the conformal weight has the form 

h^ = - + iu, h^ = - + -, (3.6) 

2 2 2 y J 



which is in agreement with the general expression ( |1.5| ). In the case II for n = odd one compares 
with the results of the polynomial solutions of the Baxter equation (TtJ to deduce that 



the quantized q% take real values. As a consequence, Uk is a real function of its arguments. 
Applying this property to ( |3.4j ) one finds that the quantum numbers q^ are real for even k and 
pure imaginary for odd k. It is for these values of the quantun numbers that one expects to find 



the state in the spectrum of the effective QCD Hamiltonian (|1.1|) with the maximal energy [11 . 



3.1. Whit ham equations 



Let us solve the quantization conditions ( [2.24 ) and ( |3.3| ) and determine the functions The 
remarkable property of the moduli Uk is that their dependence on 5j is governed by the Whitham 
equations Glj . To derive these equations one varies the both sides of the quantization conditions 



( 2.24 ) with respect to 5k and observes the following property of the differential ( |2.16| ) 



ddSo dx Ouj dx dA(x) dx 



lo d5 k y d5 k y 



N 

E 

3=3 



X 



N-j 



d5 k 



where A(x) = 2x N — x N 2 + u 3 x N 3 + ... + un and the relation ( |2.17| ) was used. This leads the 

(3.7) 



following system of Whitham equations |T4| 

du i 



05 k 



U k j(u 3} ...,u N ) 



with j = 3, N, k = 1, N — 1 and Ukj being the matrix inverse to 



v Jjk 2tt 



dxx N 3 



A 2 (x) - 4x 2N 



The Whitham equations have the following interpretation [21] in terms of "classical" reggeon 
trajectories defined by the action function Sq(x). The moduli Uk enter as parameters into the 
classical soliton wave solutions of the Hamilton- Jacobi equations. Equations ( |3.7| ) describe the 
flow of the moduli Uk in the slow "times" 5\, 5n-i, corresponding to adiabatic deformations 
of the soliton waves propagating in the system of iV reggeons. 

Analysis of the Whitham equations depends on the number of reggeons N inside the com- 
pound state. In what follows we restrict our consideration to N = 3 reggeon states and its 
generalization to higher states, iV > 4, becomes straightforward. In this case, ( |3.7| ) takes the 
form 

°' h 1 / dX , *i + S 2 = 1 , (3.8) 

'a k y(u — x)(4x 3 — x + u) 



du 2tt 
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where one integrates along two a— cycles of the elliptic curve T 3 defined in ( |2.13| ) 



T 3 : y 2 = (u-x)(Ax 3 -x + u) = -4 JJ(x - <r fe ) . (3.9) 

i=i 

The relations Q3.8| ) define as smooth functions on the complex u— plane with the only possible 
singularities at the exceptional points u = u CT a, at which two branching points merge, <jj = a^, 
and the integrand in ( |3.8p and ( |3.3| ) develops a pole. It is easy to see that these points are located 
on the real u— axis at 

11 . 
U cnt = t=, 0, — =. (3.10) 



'27 V27 

The integration in ( |3.8|) goes along a— cycles that encircle two cuts running between the branching 
points (Tj on the complex x— plane. Making the choice of the cuts one requires that the moduli 
u has to be a real function of the flow parameter S. 



3.2. Real moduli 

In order to calculate the quantum numbers q 2 and q 3 using ( |3.4j ) and ( |3.5|) it is sufficient to 
analyze ( |3.8|) only for real values of the moduli u and the flow parameters 8k- Let us specify the 
location of the branching points of the curve r 3 with real u on the complex x— plane. According 
to ( |3.9|) , three branching points satisfy the equation 4cr| — cr,- + u = and the last one is cr 4 = u. 
Introducing two parameters 

O^Oi g= 1 )2 = 1 (3n) 

0"3 -<ri' 9 V ; 12(m 2 -m+l)' V ; 



one can parameterize the branching points in the standard way as f22 

a 1 = -A(m + l), a 2 = A(2m-l), a 3 = A(2 - m) . (3.12) 
The remaining branching point, cr 4 = u, gives the value of the moduli 

1 (2-m)(2m-l)(m + l) . Q1D , 

m = -4aia 2 a 3 = j-^ . 3.13 

6V3 (m 2 -m+ l) 3 / 2 

One notes that the expression for u 2 is invariant under modular transformations, m — > 1/m and 
m — > 1 — m, induced by permutations of Oy 
Defining the effective discriminator 



A cS (u) = 16cr| Yl^i - Ojf = u 2 (l - 27m 2 



i>j 



one has to distinguish three different cases corresponding to A c fj > 0, A c g = and A e fj < 
0. Additionally using the symmetry of A e fr(w) under u — > —u one restricts the moduli to be 
nonnegative, u > 0. 

In the first case, A c g > 0, or < u < one finds that all four branching points are real 
and they can be ordered on the real axis as <j\ < < u < a 2 < cr 3 . In parameterization ( |3.11| ) 



and (|3.12|) , this corresponds to | < m < 1. The cuts of the function y(x) defined in (|3T9|) run 
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along the segments [cri,u] and [02,03] on the real axis and two a— cycles encircle them on the 
complex x— plane as shown in Fig. la. 

In the second case, A e g = 0, two branching points merge: either at 02 = 04 = for m = 1/2, 
or at cr 2 = °"3 = 7^75 for m — 1, generating two singularities on the moduli space, u = and 
u = respectively. 

In the third case, A c s < 0, or -7= < u < 00, two branching points, <J\ and 04 = u, are real 
and two remaining ones are complex conjugated to each other, o* 2 = 03. This corresponds to 
choosing the parameters as 



m 



v + i 



v — 1 



y/3i 



02 



\/3 v + 3i 



v- 



6 v^2^3 



0-3 



0o 



(3.14) 



with z/ > y/3 and the expression for moduli looks like 

1 z/(9 + i/ 2 ) 



u 



27 iy 2 - 3) 3 / 2 



(3.15) 



The restriction v > v3 comes from the condition for u to be real. Choosing the cuts between 
the branching points we require the r.h.s. of (|3.8| ) and ( |3.3D to be real. As a result, one cut runs 
along the segment [01, u] on the real axis and the second one connects the points 02 and 03 on 
the complex plane. The definition of the corresponding a— cycles is shown in Fig. lb. 







(a) (b) 

Figure 1: The definition of the a— cycles for real positive moduli: (a) < u < h== and (b) 



u > 



i_ 

27' 



In the above consideration the moduli were assumed to be nonnegative. To describe negative 
values of moduli, u < 0, we use invariance of the curve ( |3.9| ) under x — * — x and u — > —u. This 
transformation interchanges the a— cycles on the complex x— plane and leads to the relation 
u(5\,52) = —11(62,61) with 6\ + 62 = 1. Therefore, as a function of 6 = 6\ the moduli u(6) 
11 : ') . I 6) satisfies the condition 

u (l - 5) = -u(6) , (3.16) 
which maps positive and negative values of u. 



3.3. Reference points 



Before integrating the Whitham equations (3^) let us consider the quantization conditions ( |3.3|) 
at three reference points on the moduli space, u = 0, u = and u — > 00, and find the 
corresponding values of the flow parameter 6. 



15 



For u = 0, or equivalently m = |, we get from 

5(0) = 1 -I dx = 2 -f / dx =1, (3.17) 



The same result can be also deduced from (|3.16|) . The value of moduli, u(^) = 0, corresponds 



to one of the singularities ( |3.10|) of the curve T 3 , at which two cuts [<Ji,u] and [02,03] merge as 
o 2 — > u (see Fig. la). 

For u = H=, or equivalently m — 1, one approaches another singularity ( |3.10D , at which the 



cycle «2 is shrinking into a point. Two branching points merge, o<i = o 3 = and the integrand 
( |3~3"D develops a pole at x = It is compensated however by vanishing numerator leading to 



If o?x 1 f V27 dx 



(3.18) 



Applying the symmetry ( |3.16| ) one also finds that <5^— 7^ J — 0. 

For u — > oo, or equivalently z/ — > in ( |3.15| ), one rescales x — > (m/2) 1//3 x in ( |3.3|) and obtains 



the following asymptotics 



3/m^i/3J dx 6/m^V3 /- 00 dx 1/3 2vr 



iTTTF-UJ LvTTW = u ' fHlj- (3 ' 19) 

which is valid up to subleading in u terms. 

Relations ( |3.17| ), ( |3.18| ) and (|3.19| ) determine the values of the flow parameter corresponding 



to two singularities ( |3.10| ) on the moduli space and fix its asymptotics at infinity: 



u(§)=0, «(l) = -^» U{6) J^ 63 ■ (3 ' 20) 



Let us now apply the Whitham equations fl3.8|) to reconstruct the flow u = u(5) between these 
reference points. We shall consider separately two branches, < u < for | < 5 < 1 and 
u > for 5 > 1, and then glue them together at u = 



3.4. Whitham flow for < u < 4= 

V 27 



For real positive moduli inside the interval < u < h= the definition of the a— cycles is shown 



on Fig. la. The Whitham equation (|3.8| ) looks like 



85 1 f u dr 

1 (3.21) 



du n J ai y/(u - x )(4x 3 -x + u) ' 

where 5 = 5\ and the branching point oi takes real values defined in ( |3.12| ). The elliptic integral 
entering the r.h.s. of (|3.21| ) is explicitly positive thus implying that u is an increasing single 



valued function of 5. It can be expressed in terms of the elliptic function of the 1st kind defined 
in the Appendix, ( |A.1| ), as 

05 K(z) 
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Here the following notations were introduced 



(2 - m) 3 m _(m + l) 3 (m-l) 

Z ~ 7 7v\s ' (3.23) 



12(m 2 — m + l) 2 m(m — 2)' : 



the parameter m was defined in ( |3.11| ) and for the moduli inside the interval < u < -4= its 



value can be restricted as | < m < 1. 



Relation ( |3.22|) defines the derivative £ as a function of the parameter m. Replacing u 
by its expression ( |3.13|) one integrates ( p.22|) to obtain the dependence of the moduli on the 



flow parameter through the parametric dependence u = u(m) and 5 = 5(m) on the interval 
~ < m < 1. In this way, one reconstructs the branch of the curve u = u {5) for < u < and 
~ < 5 < 1 shown in Fig. 2. However, trying to obtain the explicit form of the function u = u(S) 
by inverting the relations u = u(m) and 5 = 5(m) one realizes that the parameter m is not a well 
defined function of the moduli due to invariance of ( |3.13| ) under modular transformations of m. 
Additionally, the end points, m = | and m — 1, correspond to the singularities on the moduli 
space (|3.10|) , u = and u = respectively, and one should expect to find a nonanalyticity in 
the dependence u = u(5) at their vicinity. 

In order to define a good parameter of the expansion of moduli which will replace m and 
will allow to study analyticity properties of the function u = u(5) one introduces the dual flow 
parameter 5e>. In analogy with ( |3.3| ) it is defined as an integral of the Whitham differential dSo 
on the elliptic curve T 3 along the (3— cycle encircling the interval [u, 02] 

1 p (3u-2x)dx 
Od = - 



7T J u - U)(AX Z - X + U) 

The dependence of the moduli on is described by the Whitham equation similar to ( |3.21| 

dS D 1 r 2 dx K(l-z) 



du tt J u ^(x-u)(4:X 3 -x + u) t^Vc 



(3.24) 



with z and c given by ( p. 23 ). The flow parameters 5 and 5d are not independent and their mutual 



derivative defines the Jacobi r- and q— parameters of the elliptic curve T 3 

T-i™Z- i K{1 - Z) (3 25) 

q ~ 6 ' T " 1 85 ~ 1 K(z) ' [ } 

In agreement with the Riemann theorem, r has a positive definite imaginary part and additionally 
Rer = for real moduli inside the interval < u < This property leads to < q < 1 and it 
suggests us to identify the q— parameter as a good parameter of the expansion of the solutions, 
u = u(q) and 5 = S(q), of the Baxter equation. 

3.4.1. Flow at the vicinity of u = -4= 



27 



Let us first develop the q— expansion of u at the vicinity of u = This value corresponds to 



the parameter m — 1 in ( |3.13| ) and one uses ( |3.25| ) and ( |3.23| ), combined with the properties of 



the elliptic function ( |A.5|) , to find the leading asymptotics of the elliptic parameters r and q as 



i 1 1 ~ m l-m 

t = in — - — , q = — - — as m — > 1 . (3.2o) 

tt 2 2 
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Similarly, expanding (|3.13|) , (|3.22|) and ( |3.24|) near m = 1 one gets using ( |A.2j ) 



u 



27 



'1 - m) 



85 
dm 



85 



m) 



D 



dm 



(1 — m) ln(l — m) , 

An 



where only leading terms were kept. We conclude from these relations that at the vicinity of 
u = -4= the moduli u is an analytical function of the flow parameter 5 and not of the dual 



parameter 5r>. 

It follows from ( |3.26| ) that the q— parameter vanishes at m = 1 and one can expand u and 5 
into a power series in q. Indeed, using ( |3.25|) , ( |3.23|) and ( |A.5| ) to invert the dependence q = q{m) 
near m = 1 

1 - m = 2q - 2q 2 - 2q 3 + 6g 4 - 4g 5 - 6g 6 + 0(q 7 ) 



and integrating ( 3.22 ), one obtains first few terms of the weak coupling q— expansion of the 
moduli and the flow parameter as 



u 



27 



891 



1 q z + q q 

2 H 8 H 



9 2 189 4 1917 

1 q H q 

2 H 8 16 



12879 



-q 



16 
83421 



q° + 



716499 



128 



-q 



-q 

128 

929367 
' 256 



_ 9730773 
256 

g 10 + 0(q 12 



q w + 0{q 



12\ 



(3.27) 



where only even powers of q appear. Finally, one inverts the second relation and expresses the 
moduli as a power series in the flow parameter 5 



u 



1 - 35 + 25 2 



9 81 243 V ; 



5=1-5. 



(3.28) 



This expression provides the solution to the Whitham equation (|3.21 ) at the vicinity of 5 = 1 
and u = -7=. One verifies that it coincides with analogous expression first obtained within the 



class of polynomial solutions of the Baxter equation [^, |14]]. Applying ( |3.16| ) one can also find 
the behaviour of the moduli around 5 = 0. 



3.4.2. Flow at the vicinity of u = 



The series ( |3.28| ) has a finite radius of convergence and for 5 close to 5 = ^ one expects the 
dependence u = u(5) to be drastically changed since the reference point u(|) = is another 



singularity of the curve, (|3.10| ). One finds from ( |3.13|) that u = corresponds to m = | and 



one calculates the values of the Jacobi parameters of the curve, ( |3.25| ), using ( |3.23j ) and ( |A.6| ), 
as r = and 5 = 1. Thus, to approach the point u = on the moduli space using ( |3.27| ) one has 
to develop the strong coupling expansion in q. 

Let us show that at the vicinity of u = the strong coupling expansion of moduli u in 
q = exp(i7rr) is equivalent to the weak coupling expansion in the dual parameter qo defined as 

7T 2 

q D =e~^ } \nq\nq D = —. (3.29) 

Repeating previous analysis we examine the behaviour of r and qn near m = |. One gets using 
( [OP and ( ggg ) 

in 1 2( 1\ 

to = 77 m - o ^° 



3 1n[|(m-i)]' HU 3 V 2 
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and observes that qr> vanishes at m — \ while q = 1. As m — > ~, the expansion of the Whitham 
equations (|3.22|) and ( |3.24|) using (|A.3|) leads to the following relations 



2 1 



85 
dm 



■2 In 



2 1 
- \ m 

3 V 2 



85 



D 



dm 



14 / 1 

1 m 

3 V 2 



They imply that as a function of the flow parameter u = u(8) has a singularity at 8 = \ 

as u — ► , 



A 1 3 1 

o ~ u mw , 

2 7T 



(3.30) 



and at the same time u is an analytical function of the dual flow parameter 8d- This suggests to 
search for the moduli as a function of the dual flow parameter 5d and identify the dual parameter 
q D as an appropriate parameter of the perturbative expansion of the moduli around u = 0. 
Using ( |3.29| ), ( |3.25| ), ( |3.23| ) and ( |A.6| ) to invert the dependence qo = qD(jn) near m = | 



m-\ = \q D -?>ql + \q% - Qq 7 D + - 18^ + 0{$) 



and integrating the dual Whitham equation (|3.24|), one obtains the following relations 



u 



>D 



15 

~2 
11 



459 



-9d 



6547 7 339955 □ n , 



243 



-9d - 



16 
2599 



16 



128 
113699 

128 



<& + 0(q D l 



(3.31) 



which provide a weak coupling expansion of the moduli and the dual flow parameter in qo- 
Inverting the dependence 5d = 5d(?d) we obtain from ( |3.31|) the expression for u = u(8d) as a 
power series in the dual flow parameter Sd 



u 



8 D - 28 3 D -66 5 D -488 7 D - 510 5 9 D + O (6%) 



which should be compared with (|3.28 ). Finally, to restore the dependence of the moduli on 5, 
one rewrites the definitions ( |3.25| ) and ( p.29| ) in the form 



1 3 [*> 85 D 
/ dq D \nq D -— 

2 vr J dq D 



and replaces Sd by its expression ( |3.31|) . This gives the qo~ expansion of S which can be combined 
with the first relation in ( |3.31| ) to determine the nonleading terms in the asymptotic expansion 

S ^° I _ i( u+ 2u 3 + 18m 5 ) Inu + - Uu - -u 3 - ———u^\ + 0(u 7 ) . (3.32) 

2 7T 7T \ 2 20 / 

This relation can be inverted and u can be expressed in terms of the Lambert function of 8. 



3.5. Whitham flow for u > -4= 

V 27 

The Whitham equation for u > ^= has the same form as ( |3.21| ), but important difference with 
the previous case is that two branching points, 02 and 03, take complex values. According to our 



19 



choice of the cuts between the branching points, shown in Fig. lb, the elliptic integral entering 
the r.h.s. of fl3.21|) takes real positive values defining u to be an increasing single valued function 
of 5 for 5 > 1. Since u = u(5) has the same property for < 8 < 1 it can be now extended using 
( 13. 16 ) to arbitrary real 5. 

Using parameterization of the branching points ( |3.14j ) the elliptic integral in ( |3.21| ) can be 
calculated as 



85 
du 



2^3 



v 2 -3 



71 (9 + Z/ 2 ) 3 /4(l + z/ 2)lA 



-K 



Re 



[v + i){v-3if 
_{u - i)\v + 3i) 3 _ 



1/2" 



(3.33) 



where v > \^3. Replacing u by its expression (|3.15|) one integrates numerically (|3.33| ) to restore 
the v— dependence of 5 and then calculates the moduli through the parametric dependence u = 
u{v). This gives the second branch of the curve u = u(5) for 5 > 1 shown in Fig. 2. 

To specify the appropriate boundary conditions for ( |3.33| ) one considers two values of the 
parameter: v = \/3 and v — > oo. For v — > \/3 one finds from (|3.15| ) that u infinitely increases 
and its leading asymptotic behaviour in 5 is given by ( p.20|) . For v —>■ oo one gets from (|3.15|) 
and ( |3.20| ) the corresponding value of moduli as«= ^= for 5 = 1. Let us consider the flow of u 
in both cases in more detail. 

For large v one uses ( |3.15| ) and ( |3.33| ) to expand u and 5 in inverse powers of 1/V 
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Inverting these relations one obtains the expansion of u in powers of 5 — 1 — 5, which identically 
coincides with (|3.28|) . This means that two different branches of the function u(5) corresponding 
to the Whitham flow for | < 8 < 1 and 5 > 1 can be smoothly glued at 5 = 1 (see Fig. 2). 
For v — > \/3 the moduli ( p . 1 5[ ) goes to infinity as 



u 



A(u 2 - 3) 



-3/2 + {u 2 



3)- 1 / 2 + 0{{v 2 -3f 2 ). 



To find the corresponding behaviour of the flow parameter one takes into account that S 
v = oo and integrates the Whitham equation ( |3.33|) as 



(3.34) 
= 1 for 



5 = 1 + 



v 1 du 85 
dv 

oo 



dv du 



Replacing the derivative by its expression Q3.33Q , we expand the integral around v 2 = 3 and 
apply (|A.4|) to get after some calculations 

2 5 / 3 7T , , 



5 



r 3 (l) 



3) 



-1/2 



l + 0{{u 2 -3f/ 2 ). 



(3.35) 



Combining together ( |3.34| ) and (|3.35| ) we obtain the asymptotic behaviour of the moduli at large 
positive values of the flow parameter 5 
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r 9 (|) 

(2tt) 3 



(5 3 + 35 2 ) x [1 + 0{l/5)\ , as 5 -> oo . 



(3.36) 
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Here, the leading term coincides with (|3.20|) . The symmetry ( |3.16|) allows to extend the flow to 
large negative u as 



u 



r 9 (l) 

(2tt) 3 



(5 3 - Q5 2 ) x [1 + 0(1/ 6)] , as5^-oo. 



(3.37) 




Figure 2: Quantized moduli u for real values of the flow parameter 5. The dots indicate the 
positions of three singularities of the elliptic curve 1^ located at 5 — \ = ±| and 0. The 
behaviour of the function u = u(5) around these points and at infinity is described by (|3.32|) , 
( CT and (£3 



Summarizing the results of this Section we present on Fig. 2 the dependence of the quantized 
moduli u = u(5) on the flow parameter — oo < 5 < oo. According to (|3.16| ) the moduli u is an 
odd function of 5 — |. Its asymptotics around S = 1/2, 1 and oo is described by ( |3.32| ), (|3.28|) and 
( |3.36|) , respectively. We recall that this curve provides the solution to the quantization conditions 
rSf) and (|3.S|) for iV = 3 reggeon compound states. 



4. Quantum numbers of N = 3 reggeon states 

Let us apply the results of the previous section to evaluate the quantum numbers q 2 and q 3 of 
the N = 3 reggeon states in the leading order of the WKB expansion. According to (|2.25 ), the 
values of qs depend on two complex numbers £\ — ^ -\-iv\ and £2 = ^ + iZy 2- Following ( |3.4| ) and 
( p.5|) we consider two special cases, n\ = n 2 = and u% — i/ 2 — 0, corresponding to real values 
of the moduli and the flow parameters.^ In the case I, the relations ( |3.4D and ( p.6[ ) between the 
moduli and quantized q 3 can be written in two equivalent forms 



\3 

[IV) x u 
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or 
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-1 



^1 



Vi + v 2 



(4.1) 



(4.2) 



lr To obtain real values of u and 8 for iV = 3 states it is enough to impose even weaker condition: ni/712 = v\/v2- 
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with ui t 2 being real. In the case II, one gets similar relations by replacing iv^ — > ^n^ and taking 
ni >2 to be integer. 

It is easy to see from ( |3.16| ) that q$ is antisymmetric under interchanging of V\ and z/ 2 , while 
the conformal weight h = | + i{v\ + u 2 ) is explicitly symmetric. Therefore, q% vanishes for 
Vi = z/ 2 and the corresponding N = 3 reggeon state with the conformal weight h — | + 2iv\ is 
degenerate f9|. Their wave function does not depend on one of the reggeon coordinates and the 
corresponding energy is equal to the energy of iV = 2 reggeon state, E 3 (q 3 = 0, h) — E 2 {h). 

Since quantized q 3 and h depend on the same parameters v\ and v 2 one can express z/ 2 in 
terms of h and consider q 3 to be a function of the conformal weight h and real v\ . Let us consider 
separately the dependence of q 3 on its arguments, q 3 = q 3 (h, v\). 

For fixed values of the conformal weight h = fixed it is convenient to apply ( [4 . 1| ) . One finds 
that up to rescaling of the argument the ^—dependence of q 3 is governed by the moduli and it 
follows the pattern of the curve shown in Fig. 2. 

For fixed values of v\ one uses instead Q4.2| ). Introducing notations for q = u(5) 5~ 3 and 
x = 5~ l one writes (14.21) as 



?3 



{iviY x q(x) 



h- 1 - 
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lUl X X 



q = x 3 u(l/x) 



(4.3) 



with q depending on x and not on v\. Here, the parameter v\ fixes the scale of q 3 and h—\. It 
takes real continuous values which cannot be however arbitrary small for the leading order WKB 
approximation to be applicable. In the case II one replaces iv\ — > ^rix and uses integer rii to 
label the curves ^3 = qs(h; n-i). Again, the values of n\ cannot be small. 

The function q = q(x) can be determined out of the curve u = u(8) depicted on Fig. 2 in 
two steps. One first finds the dependence of moduli u on x = 1/5 as shown in Fig. 3 and then 
obtains q = q(x) in the form shown in Fig. 4. 



. 6 T 



0.4 




-0V2 

^\ -0.4 
\ -0 . 6 



Figure 3: Two different branches of the function u(l/x) entering into (|4.3| ). Dotted lines represent 
u = ±77^ and two dots indicate the singular points x = 1, u = ^7= and x = 2, u = 0. The part 



of the curve u 2 (x) < corresponds to the polynomial solutions of the Baxter equation and it 
should be compared with Fig. 4 in [|T0 . 



The following comments are in order. To describe the transition from Fig. 2 to Fig. 3 one has 
to take into account that 5 = is mapped into infinities x = ±00 and two parts of the continuous 
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function u = u(6) for 5 < and 5 > give rise to two branches of the functions u and q on 
Figs. 3 and 4, respectively. 
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Figure 4: The function q = q(x) defining the quantum numbers ( [4.3|) of the N = 3 reggeon 
compound states. Its behaviour around x = 0±, x = 1, x = 2 and x = ±oo is given by (|4.5|), 
(p|), (p|), O and (0), respectively. 



The asymptotics of q as x — > ±oo follows from ( |3.28| ) and ( |3.16j ) 
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as x — > ±oo . 



(4.4) 



Two infinities 5 — > ±oo are mapped into the origin x = 0. Since u = u(5) has different sub- 
leading terms in the asymptotics ( |3.36|) and ( |3.37| ), the function q has different asymptotics as x 
approaches the origin from different sides (see Fig. 4) 
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It is this property that is responsible for the appearance of a cusp on Fig. 4 at x = and 



q = T ^t)? = 0.0617... Two reference points ( 3.17 ) and ( 3.18 ) correspond to x = 2 and x = 1, 



respectively, and the asymptotics of q at their vicinity can be found from ( |3.30| ) and (|3.28|) as 
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(4.7) 



The part of the function u = u(5) corresponding to < 5 < 1 and u 2 < describes the 
polynomial solutions of the Baxter equation. It is mapped into the branch that starts at q = 
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for x = 1 on Fig. 4 and then decreases to infinity according to (|4.4| ) as x increases. One checks 
that this behaviour is an agreement with the numerical results (see Fig. 4 in [|i~0fl ). 

We would like to stress that the above results were obtained in the leading order of the WKB 
expansion and they are valid for large values of quantum numbers. Examining Fig. 4 and using 
( |4.3| ) one notices that the latter condition is violated at x = and x = 2 when either \h — || or 
q 3 vanish. Therefore, one should expect that the behaviour of the function q = q(x) close to the 
cusp on Fig. 4 and around the point q = will be modified by the nonleading WKB corrections. 
Their detailed analysis deserves further investigation. 




Figure 5: The energy spectrum of the N = 3 reggeon compound states evaluated in the leading 
order of the WKB quantization for type-II quantization conditions, E 3 = E 3 (h;ni). Three 
different curves correspond to n\ = 4, 8 and 12. The dotted line represents the BFKL energy of 
the N = 2 reggeon state, E 2 (h) = 4 [#(1) - *(| + \h - ||)] . 

Concluding our consideration we would like to apply the obtained results to calculate the 
energy of the iV = 3— reggeon states. The corresponding expression looks like 

3 

£ 3 WKB fe, qs) = -2\n2-J2 W + + *(1 " - 2*(1)] (4.9) 

3=1 

where = ^lnT(x) and Sj are defined as zeros of A(x) 

A(5 j ) = 26? + q 2 6 j + q 3 = 0. 

One should warn however that, first, this expression was obtained for the type-II solutions (|3.5| ) 
and, second, it is not exact and is valid in the leading order of the WKB expansion. Replacing 
iv\ — > ni/2 in (f4.3|) and substituting quantized q 2 and h into (|4.9|) one gets E 3 = E^{h;ni) as 



a function of the conformal weight h and integer ri\. To get an insight into the spectrum of 
i?3 we plot on Fig. 5 the dependence of E% on the conformal weight h for different values of 
the integer n\. These curves are in an agreement with the numerical results of the polynomial 
solutions (see Fig. 8 in 0]) and they support the calculation of the odderon intercept performed 
IIJ. For given n\ the energy is maximal for h = ^ and its value increases as n\ decreases. The 



in 



absolute maximum of the energy Eq = max^^E^/i; n\) defines the odderon intercept and it is 
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expecting to occur at ri\ — 1 and h — |. One also observes that for the upper curve on Fig. 5 
the maximal energy, E^{\\ 4) < i?o, is close to the lower bound on the odderon energy obtained 
in the variational approach pi]. We refer to |K| for detailed description of the properties of 
curves shown on Fig. 5. It would be interesting to improve (O) by applying the expression for 



the energy E 3 proposed in [11 



5. Summary 

In this paper, we have studied the spectrum of the reggeon compound states in high-energy QCD. 
These states appear as solutions of the (2 + 1)— dimensional Schrodinger equation ( |1 . 1|) which 
exhibits remarkable properties of the SL(2,€) invariance and complete integrability. As a result, 
the energy of the states can be evaluated as a function of the set of quantum numbers q 2 , qN- 
The latter are defined as eigenvalues of the mutually commuting integrals of motions and their 
possible values are constrained by the quantization conditions. 

We have established the quantization conditions by applying the methods of nonlinear WKB 
analysis to the N— reggeon Schrodinger equation in the separated coordinates. In the leading 
order of the WKB expansion, the N reggeon state looks like the classical system of N particles 
moving on the 2— dimensional plane of transverse gluon coordinates along periodic trajectories. 
This collective motion can be also considered as a propagation of the soliton wave in the system 
of N reggeons. The charges q 2 , q^ become parameters of the soliton waves. 

The quantization conditions follow from the analysis of the wave function of the N reggeon 
state in the leading order of the WKB expansion. Requiring the wave function to be a single- 
valued function of the separated coordinates, we have found that the selection rules for q 2 , 

qN have the form of the Bohr-Sommerfeld quantization conditions imposed on the reggeon 
classical trajectories in the action-angle variables. The same conditions can be interpreted as 
Whitham equations on the moduli of the spectral curve r^r corresponding to the classical N 
reggeon system. We have solved the Whitham equation for N = 3 reggeon states by using the 
properties of the elliptic integrals on the curve T 3 . The quantized q 2 and q 3 were obtained in 
the form of perturbative expansion in powers of the Jacobi q— parameter of the elliptic curve 
Different parts of the spectrum of q 2 and q% correspond to the weak, g C 1, and strong, 
q ~ 1, coupling regime of the perturbative q— series. In the latter case, one performs the duality 
transformation to express q 2 and q% as a weak coupling expansion in the dual coupling constant, 
qn <C 1. Combining together different branches we have obtained the spectrum of quantized q 2 
and #3 which is in agreement with available exact solutions found within algebraic Bethe Ansatz 
approach. 

In conclusion, one should mention that the above consideration was restricted to the leading 
order WKB expansion and one should additionally study the importance of nonleading correc- 
tions. In particular, considering the behaviour of ( |2.10| ) around the origin x = one can argue 
P~0f| that nonleading WKB corrections to Q(x) become equally important. As in the case of the 
polynomial solutions, ( j2.9p , the analysis of the singularities of the solutions to the Baxter equa- 
tion at x = should lead to additional constaints on the quantum numbers of the reggeon states 
and presumably fix the ratio C + _/C_ + in (|2.8|). These questions deserve additional studies. 
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Appendix: Elliptic integral of the first kind 

In this appendix we collect some useful properties of the elliptic integral of the first kind. It is 
defined as |22], p3| 

K(z) = C dX = = -F(ii; 1; z) (A.l) 

1 ' Jo v /(l-a;2)(l_^2) 2 ^ 2 ' 2 ' ' y 1 

with F being a hypergeometric function. For z — > it behaves as 

K « = i{ 1 + G) , * + £i)'* , + -> (A2) 

and for z — > 1 as 

K(z) = -iK(l - *) In 1^ - 2(|) 2 i(l - z) - 2(M) 2 (^ + s^)(l - + - (A.3) 
Calculating ( |3.35| ) one uses the relation 
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The Jacobi r and q parameters defined in ( |3.25| ) have the following asymptotics for z — > 

i. 16 i / 13 2 23 3 2701 4 /n/ 5 , 

r(z) = -In [z-\ z 2 H z H z A + 0(z 5 

K ' 7T z 2tt V 32 96 16384 v ' 

x fi 1 21 2 31 3 6257 4 /n/ 5 ,\ 

= — l + -zH ,2 H z 3 H z 4 + £>(z 5 ) . (A.5) 

qK 1 16 V 2 64 128 32768 / 

and their asymptotics for z — > 1 can be found using the relations 

r(z) = - ./ r , = exp f — ^ -V (A.6) 



r(l — z) ' \ln g(l — z 

which follow from the definition (|3.25|) . 
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